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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


INTRODUCTION 


In a previous article [1] the authors defined in geometric terms 
the decoupling problem for a constant linear multivariable system: 
namely, the problem of achieving independent control of specified 
outputs by means of suitably combined inputs and of suitable linear 
state variable feedback. Necessary and sufficient conditions for 
decoupling to be possible were found in two important cases; but the 1 
general problem is unsolved. However, if in addition to state feed- 
back, dynamic (integrator) compensation may be utilized, it becomes 
possible to state general necessary and sufficient conditions for ^ 

decoupling in a simple and constructive way. Geometrically the decoupl- 
ing synthesis amounts to extending the state space of the original 
system to a larger space, the increase in dimension being the number 
of integrators used in dynamic compensation. In addition, state- 
space extension can be used to achieve a desired pole distribution for 
the closed loop system transfer matrix. 

In the present article we state and solve the extended decoupling 
problem (§ 1). Under certain restrictions, the problem of minimiz- 

1 

ing the order of dynamic compensation (i.e., the dimension of the ex- 
tended state space) is solved in § 2 . This solution is actually , 

the best possible if the number of scalar inputs is equal to the 
number of output blocks to be decoupled (§3). In § 4 the role 
of state space extension in pole assignment is determined. It is 
shown that, with dynamic compensation of high enough order, any pole 
distribution can be synthesized for the decoupled system, whenever 
decoupling is possible at all. An example is given in § 5. In 
conclusion (§ 6) a more general view of decoupling is taken, with 
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the restriction to linear compensation relaxed. The resulting open 
loop decoupling problem is shown to be equivalent, however, to the 
extended decoupling problem of § 1. 

In the sequel, the material in [1] is assumed to be known. 

NOTATION 


Script letters E,E',R,N,... denote vector spaces over the reals, 
with elements x,y,...; d(E) is the dimension of E; U ■ 1/ means U, 1/ are 
iomorphic, i.e., d(U) = d(l/). A,B,C,... are linear maps; A|R is the 
restriction of A to R; B or {B} is the range of B. Spectrum means com- 
plex spectrum. A symmetric set of complex numbers is one of the form 

^1 1 ' 8 1 1 6 1 » 02 » $2 ' •••} 

where the ou are real and 8^ is the complex conjugate of 0^ . W (H) is 
the kernel (null space) of H. 

With A,B,E' fixed, C(l/) is the set of maps C such that (A+BC) V 
CV, C' (V) the set of C such that (A+ (B+E' )C) 1/ c l/. I_ (resp. ) 
is the class of V such that C(l/)^<)> (resp. C' (1/)^<J>). If d(E) = n, 

A : E-*E and BCE, then 



(A | B} 



RCE is a controllability subspace (c.s.) for the pair (A,B) , 
written Re C, if C(R)^4> and if, for some Ce C(R), 


R = {A+BC | BOR} ; 


R is determined uniquely, as written, by any Ce C(R). Similarly S is 


2 - 



a c.s. fcr (A , B+E 1 ) , written Se £' , if C' (S)^4> and 
S = {A+ (B+E ' ) C | (B+E’)nS}, Ce C'(S) 


The maximal (i.e., largest) element of I_ (resp. £) contained in a 
subspace T is denoted by max (I_,T) (resp. max (£,T)), and similarly 
for V , £' . It is known from [1] that these maximal elements exist 
and are unique for each fixed T and that, if 1/ = max (I ,T) , then 

max (£,T) = {A + BC | BOl/}, Ce C(l/) 


J is the set of integers (l,...,k). Unless otherwise noted, 
all summations and intersections are over J. If R^, ieJ, is a 
family of subspaces, 


R* = HR., 

l -Hi D 




R* = n r. 
i 1 


A[R., J] 5 J d(R i ) - d(E Ri> 


Certain auxiliary results needed are collected in the Appendix, 


1. EXTENDED DECOUPLING PROBLEM 


As in [1] the control system is specified by the differential 
equation 

x (t ) = Ax (t ) + Bu (t) (1.1) 

and output relations 

y i (t) = H ± x (t) ieJ (1.2) 

The state vector xeE, d(E) = n ; the control vector ueli, d(U) = m ; 
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the output vector d(y^) = q^. The maps A , B , H ^ are independent 

of t; in fact, (1.1) qua differential equation plays no role until 
§ 6 , as our problem is purely algebraic. 

Write = W(H^), ieJ; as in [1] we assume AT f E, ieJ. In [1] 
we discussed the restricted decoupling problem (RDP) : Find R^e£ (ieJ) 

such that 

0 C(R i )^4> (1.3) 

R. c n W ieJ (1,4) 

1 37*i 3 

R i + W i = E ieJ ' (1.5) 

A family of c.s. R^ e £ (ieJ) which satisfies (1.4) (but not necessarily 

(1.3) or (1.5)) is admissible . Let R^ be the maximal admissible c.s. 

M 

was denoted by R^ in [1]). It is clear that RDP is solvable only 
if 

R i + W i = E ieJ (1.6) 

In general (1.6) is not sufficient for solvability of RDP because (1.3) 
may fail for the R^, i.e., there may not exist any C such that 

(A + BC)R^ c R^ , ieJ. To avoid this difficulty we introduce an 
.ctended decoupling problem as follows. 

Adjoin to (1.1) the equation of a new dynamic element: 

x' (t) = I'u' (t) (1.7) 

where x'eE', u'eU\ d(E') = d(U') = n' , and I ':U'zE' ; 

the input u'(*) can be freely chosen. For the system (1.1) extended 
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by (1.7) define the state space 


E e = E © E ' 


and the extended input space 


U e = U © IT 


Define extensions A e , B e , E' of A,B,I' as follows 


A e : E e ->-E e ; A e (x + x') = Ax (xeE, x'eE’) 


B e : U e +E e ; B e (u + u') = Bu (uetl, u'eli') 


E': u -*-E e ; E' (u + u') = I'u'fueU, u'eli') 


( 1 . 8 ) 


6 6 6 

Below we write A, B for A ,B ; x for vectors in E ; and P for the 
projection E©E'-*-E. Observe that PA = AP = A, PB = B , PE* = 0. 

The combined system (1.1), (1.7) is now specified by the pair (A,B+E') 
The extended decoupling problem (EDP) is the following: Given 
the original maps A :E-*-E, B:(i-*-E, and E (ieJ) , find:- (i) E' ( that is , 
n'), (ii) extensions A,B,E' as in (1.8) , (iii) e C' (ieJ) , with the 
properties 


n C' (S.)*4> 


s i c 


(1.9) 


( 1.10 


Si + (W ± © E ' ) = E©E' 


( 1.11 


It is clear that the choice of isomorphism I', and so of E' in (1.8), 
can be arbitrary after n' is fixed: for instance, I' = n'xn' 
identity matrix, in the coordinates selected. 



EDP has the same structure in E e as RDP has in E, but flexibility 
is gained from the special form of the new system map A and constraint 
spaces © E ' . Justification of EDP as the correct description of 
decoupling by dynamic compensation is clear: the output relations (1.2) 
are preserved on replacing by N^©E’ (equivalently by defining 
extensions FL of to be zero on E 1 ) ; no additional control inputs 
(B) to the original system (1.1) are postulated; subject to the latter 
constraint, full linear coupling is allowed between (1.1) and (1.7). 

Our main result (Th. 1.1) states that decoupling by dynamic com- 
pensation is possible if and only if the maximal admissible c.s. R^ 1 
of RDP are sufficiently large. 

Theorem 1.1 

M 

For the RDP of (1.3) - (1.5) , let R^ be the maximal admissible 
c.s. in £. The corresponding EDP of (1.9) - (1.11) is solvable if and 
only if 

R i + W i = F ieJ (1.6 bis) 

Proof 

1. (Only if) We show first that S e £’ implies R = PS e £. 

Since £' (S)?*<t>, ASCS+B+E', and AR = PASCR+B, so that C(R)?<p. 

Also, by Th. 2.1 of [1], S = lim S y (y = 0,1,2,...) where S° = 0 , 
o y+1 = S n (AS y +8+E ' ) . Write R y = PS y . Since PE’ = 0, Prop. 

A. 4 implies 

R y+! = p5 P+l = RH(AR U + B) ; 

again by [1], Th. 2.1, R = lim R y e£. Thus S^e£* (ieJ) implies 
PS i eC (ieJ) ; and (1.10), (1.11) yield 
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ieJ 


( 1 . 12 ) 


PS i C w -i 

1 1)7*1 3 


PS i + w i = E 


ieJ 


(1.13) 


By (1.12) and maximality of the , PS^^C R^ ; this and (1.13) imply 

( 1 . 6 ). 

2. (If) Assuming (1.6) holds, define n' = £d(R. M ). With n’ so 

i 1 

large, there clearly exist maps M^:E e -► E ' with the properties: 


R^ONfM.) = 0 , {M i > = M i R^ 


ieJ 


and the ranges {M^} (ieJ) independent. Define = (P + M^)R^ M (itJ) . 
Then 


AS. = AR“ c + B C 5 i + B + E', ieJ ? 

t 

and since the are clearly independent, there exists C: E e -*-(/ e with 
Ce 0 C(Sj)» It will be shown that S^e£’ . Dropping the subscript i, 
suppose ReC, so that the relations 

R° = 0, R m+1 = RO(AR u + B) (y = 0,1,..) 

imply R^TR. Let tM}CE» and 

S = (P + M)R, S° = 0, S y+1 = SO(AS U + 8 + E') 

S° (P + M) R° ; and if S U D (P + M) R y , 

s y+1 D [ (P + M)R] n [A(p + M)R U + B +E ' ] 

= [ (P + M)R] n [AR U + B + E’ ] 


Then 


D (P + M)[Rfl (a;< U + 6 - E' ) ] 


= (P + M)R U+1 

By induction S D S p D (P + M) R M T (P + M) R = S, i.e. , S U fS ; so SrC' . 
Application of this argument to the R^ and yields the desired 
result . 

The relation PS. = implies 

1 .1 



n (N. 

3 


© f) 


(1.10 bis) 


By (1.6) 


S t + (I + M i )W i = (I + M i )E 

and addition of E' to both sides yields (1.11). $ 

Remark 1 The proof reveals the symmetry between a c.s. and 

its extension: if ReC and S ? (I + M) R with {M}CE' then Se£' . 

Conversely if Se£’ then R = PSe£. 

Remark 2 In part 2 of the proof the S^ were constructed to 

be independent. By [1], Th. 2.2, Cc n C(S ) can be chosen such 

i “ i 

that, for each i, the spectrum of (A + (B + E')C) jS^ is any symmetric 
set of d(S^) complex numbers. 

Remark 3 Condition (1.6) is not implied by controllability of 
/ B) , i.e., by the condition { A | B } = E. For example let 



0 10 


k* m 

0 0 

A = 

0 0 1 

, B = 

1 0 


0 0 0 


0 1 




L 

H i - 

(1,0,0) , 

H 2 = (0 ; 1,0) 
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By the methods of [1] one finds 



and (1.6) fails for i = 1, 2. 

The following description of the structure of a decoupled system 

will be applied later to solutions of EDP. The result is stated for 

RDP for simplicity of notation. Let R^, ieJ, be any solution of RDP , 

write R = iL R . , and write A = A + BC for CeC i Oc(R.). Let x be 
^ l c — 1—1 

the coset of x in E/R*. Noting that A c R*CR*, we define the induced 

map A : E/R*-*>E/R* by A x = A x. 
c c c 

Theorem 1 . 2 

A 

There exist R.cR.» ieJ, independent of CeC such that 

i ^ i ' * 

R = R*© R, © ... 0 R k (1.15) 

and 

a 

R ± s TSf. = (R i * R*)/R* (1.16) 

A 

The satisfy 

= E ieJ (1.17) 

AR.cR.©R* ieJ, CeC (1.18) 

T he spectrum of (ieJ) can be assigned as any symmetric set of 

d(!T) complex numbers by suitable choice of CeC. 

>N A 

Proof Let R^ be any subspace such that R^ = R. ® R.^R*. Independence 

A 

of R*, R^ (ieJ) follows by Prop. A.l; hence (1.15) is true, and (1.16) 
is clear. Since 
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R* C 


n 


n L 

i j^i a^j'a 


.W = 0 W. ( 


(1.17) follows from (1.4) and (1.5). 

Let C eC be fixed; write a = A 
o — o 

projection: E+E/R* ; thus A q Q = QA q . 

range BHR^, and write = QB^, IT = 

E. is a c.s. for the pair (A , B . ) . 

1 O 1 


Since A R. c R. , (1.18) is clear 
ci i 

_ , A = A ? and let Q be the 
c ° c 
o o 

Let B^:U-»-E be any map with 
Q(B0R i ) . It will be shown that 
In fact 


T. = QR. = Q(A | BOR.} 

1 1 o 1 

= {A I Q(BnRJ) 
o 1 1 

■ a o I V . 

and the assertion follows. By Prop. A.l the are independent. 

Hence (cf. [1] §§ 4,5) there exist : E / R*-*-U (ieJ) such that 

Dl. = 0 (i, jeJ ; j^i) , (S + B.D.)R.Cff. (ieJ) and (A + B.D. )|R. 
ij o iii i o ii'i 

(ieJ) has any pre-assigned spectrum. Define - D^Q (isJ) . Then 
(R j + R*) = 0 (i,jeJ ; j^i) and B.DXc B.cR.. Let D:E-*-U be any 
map such that BD = S B.D. ; D exists since (B.}CB (ieJ) . Then the 

iii l 

map C:E-*-U defined by C = C q + D has the properties required. 

2. MINIMAL STATE SPACE EXTENSION 



Theorem 1.1 shows that if (1.6) holds, EDP can always be solved 

by dynamic compensation of order n' < d(R"). There is then a 

least integer n >0 for which EDP is solvable with n' = n : in case 
o — o 

n Q = 0, the corresponding EDP reduces to RDP . From a practical view- 
point it is of interest to find n Q : we call this the problem of 

minimal state space extension, or of minimal solution of EDP. 
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The general problem of minimal extension includes the general 
solvability problem for RDP, and is unsolved. However, suppose the 
additional constraint is imposed, that 

P S L = R*? iej (2.1) 

M 

where the R^ are the maximal admissible c.s. in £. In this case it 

will be shown how to compute the minimal n' , say n M « In general 

n M >n Q , because (2.1) rules out extension of any R^t£ which is properly 
M 

contained in R^ , but which still may be large enough to satisfy (1.5). 
However, if d(8) = k, it will be shown in § 3 that (2.1) holds for 
any solution of EDP, hence n w = n , and so this case will be solved 
completely. 

It is convenient for later purposes to adjoin to (2.1) the 
additional constraint 

PS* Cl/ (2.2) 

where 1/ is a subspace such that 

l/e£, l/c(R M )* (2.3) 

In (2.3) (R M )* is the * space (see Notation) of the family R^ , icJ. 

Relations of form (2.2) arise in the synthesis of pole distributions 
(§ 4). With 1/ fixed, let 

s n E (t"nV) (2.4) 

° i j*i 3 

M 

In the remainder of this section we write R^ = R^ , ie{0}u J. 
Theorem 2 . 1 

For the RDP of (1.3) - (1.5) let R^ (ieJ) be the maximal 
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admissible c.s. in C, and assume (1.6) is true . If V satisfies (2.3) 
and if 

d(E’)>n M (lO e A[(R i + R 0 (l/))/R 0 (l/) f J] (2.5) 

then a solution S^ (ieJ) of EDP exists such that PS^ = R^ (ieJ) and 

S* C R (l/). 
o 

Conversely if EDP has a solution S^ (ieJ) such that PS^ = R^ (ieJ) 

then 

PS*eI_ , PS* CR* (2.6) 

If for some \J , PS* C V then PS* C R q ((/) and (2.5) is true. If equality 
holds in (2.5) then PS* = R q ( 1/) and (S^ + S*)nE' = 0 , ieJ. 

Corollarv 

For the RDP of (1.3) - (1.5) , suppose (1.6) is true and let 
M 

l ' = max (I_, R*) . Under the constraint (2.1) there exists a solution 

{F',S^, ieJ} of EDP if and only if d (E r ) > n M (U M ) . 

Existence of S^ will be proved by a refinement of the construc- 
tion used in the proof of Th. 1.1. For this we need Lemmas 2.1 - 2.3. 

Of these the first two assert general properties of extensions. 

Lemma 2 . 1 

Let ILC E (ieJ) . rf d(E')^ 6 = A[IL, j] there exist maps 
• E e -*E ' (ieJ) such that the subspaces l/^ = (P + M^)IK (ieJ) are 
independent . 

i-1 

Proof Write W. = 0, W. = U . n E U . (i = 2, . . . ,k) . Then 
J- ii j=l 3 

k k i-1 

dyW ) = £ Cd (U . ) + d( £ U J- 
1 i=2 1 j=l J 

~ d( £ d .) ] 


E 

i=l 
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hence there exist such that W(M.)niil, = 0, and the 

(M^ } , ieJ, are independent. Suppose the are not independent and let 

i^2 be the greatest integer such that 1/^0 ^ 7 * 0. There is x / 0 
such that 


where 



i-1 


x = (P + M i )u i 

ii 

(P + M.) 

u.eli. (l<j<i), 

so that 

i-1 


i-1 

Pu i = u i = X 

3 = 1 

PU. = 

E u j 

j=l 3 


and u^eW^. By independence of the {M^ } , M^u^ = 0, hence u^ = 0 and 
x = 0 , a contradiction. || 


Lemma 2 . 2 


Let V , R^c E (ieJ) and define 

R 5 0 S. (R. n V) 
o i 37*i j 


S = A[(R i + R o )/R Qf J] 


If d(E')>6 there exist maps NL : E e ->-E' (ieJ) such that , if 
l/ i = (P + M i )R i (ieJ) , then V* = R . 

Proof Write IT. = (R. + R )/R (ieJ) and let P be the projection: 
1 l o o 

E © E '-*■ (E/R q ) (?) E ’ . By Lemma 2.1 there exist M^: (E/R q ) © E '+E ' such 
that F^ e (P + M^)fT (ieJ) are independent subspaces of (E/R q ) ©E’. 
Let = M^P ; then IT is well defined. Since = (l/^ + R q )/R q , it 
follows by independence of the F^ (ieJ) and Prop. A.l that R q DI/*. 
For the reverse inclusion observe that, by (A.l) and (A. 3), 


and that M^R Q = 0 (ieJ) - Thenx eR Q implies x = with 

x. = 2*i x • • / x • e R . n R ; x • • e R • n R 
1 13 1 ' o A i;j j ' 1 o 

. and 

x . = (P + M.) Xi e V L ; = (P + M.)x i;j e V. 

Thus x i el/ i n l/J, so x e l/* . | 

amma 2 . 3 

Let R^ (ieJ) satisfy the hypotheses of Th, 2.1 and let 1/ satisfy 
(2.3) . If R q is defined by (2.4) then R q eI_. 


Proof Since R* c (ieJ) there follows 1/cV. (ieJ) where 

IT = max (I_, jQ^Wj). Hence for each ieJ there exists C^:U-*-E in 
C(l/^) n C (l/) . Since R^ = max (£, l/^) , there follows * 

R i = {A + BC i | 80 lM, so that C jz CfR^h C(l/)c C(R i 0 l/) . That is, 
R.nVel (ieJ) , hence £ (R.nV)el (ieJ) . Now apply the same argument 

i* 1 3 - — - 

to the pair of subspaces R^, l/^ = 2^. (Rj^l/) to get that R^n e I_ 

(ieJ) . Finally, use (A.l), (A. 3) to obtain R = £ R. ^ I/, e J_. 

O ^ 1 i- 

Proof of Theorem 2 . 1 (direct statement) Lemmas 2.2 and 2.3 provide 
and l/^C E@E* (ieJ) with the properties: d(E') = n M (U) » and 


Pl/ i = R i (ieJ) ' v * = R 0 ' ^ E i 


(2.7 a,b,c) 


Since ]_C.V , V* e I_’ . Also, by (2.7a) 


AU i c A(R i + E') c R ± + 8 C l/ i + E' + B, ieJ; 


hence 1/. e I 1 , and 
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V. + I/* e r 
1 - 


ieJ 


( 2 . 8 ) 


Because the factor spaces (l/^ + (/*)/(/* are independent, there exists 

C e H c' (V. + l/*) . Define 

i ~ ! 

S L = {A + (B + E')C | (8 + E')n(l/. + l/*)} ieJ (2.9) 

It will be shown that PS^ = R^ and S*CR q . By Remark 1 after Th. 1.1, 
PS . e C ; also 

l — 

PS i c Pd^ + V*) = R ± + R q c R. + R* C j p i N. 

Since R^ is maximal, PS^cR^. For the reverse inclusion, by Prop. A. 5, 
PS i O P[ (8 + E , )0(l/ i + (/*) ] = BO(R. + (/*) D BOR. 

Since PS^C R i there exists Ch e C(PS i )OC(R i ) . Thus 

PSj^ = (A + BC i | 8 npS i ) 3 {A + BC i | BflR.} = R. ; 

and so PS^ = R^ (ieJ) . Finally, by (A. 2), 

S* C 0 5. (I/. + l/*) = l/* = R 1 

The idea of this proof was to use (2.9) to manufacture 
'compatible' c.s. contained in the l + V*. The method works because 
the + [/* satisfy (2.8). For this one needs (2.7c), which is 
guaranteed (Lemma 2.3) by maximality of the R^, and also , which 

follows by R^ e I_. Maximality ensures also that R^DPS^. 


8 

a 
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Proof of Theorem 2.1 (converse statement) Since S^ (isJ) is a solution 
of EDP, S*e V , hence PS*eI_. Since PS^ = R^ (ieJ) , clearly PS*CR*, 
so (2.6) is true. Let PS* Cl/. Then S* = S ^ O S* (ieJ) implies 
PS*C ^ i (RjOl/) (ieJ) and so PS*cR q ( 1/). By Prop. A. 2 (where S Q is 
defined) 

d(E , )Z6 1 i A[(S. + S q + E')/(S o + E'), J] 

= A[r (S. + S )/PS„ , J] 
l o o 

= A[(R i + R q )/R o , J] 

Finally, if d(E r ) = 6^ Prop. A. 3 implies S* + E' = S q , hence 

PS* = R ; and also (S. + S*)OE' = 0. I 

o 1 ■ 

Proof of Corollary Any solution of EDP subject to (2.1) satisfies 
(2.6), hence PS*Cl/ M , and by (2.5) d ( E ' )Zn u ( V M ) . Thus n M (l/ M ) is the 
least integer for which EDP is solvable subject to (2.1). || 

3. MINIMAL EXTENSION WHEN d(8) = k 

Assume d(8) = k and let S^ (ieJ) be any solution of EDP. It 
will be shown that 

= PS i = R** ieJ (3.1) 

By Remark 1 after Th. 1.1, R^ e C and clearly the R^ satisfy (1.4), 
(1.5). It is enough to show that 

d(BnR^) = 1 ieJ (3.2) 

In fact, since ^ E, (1.5) implies R^O, hence (3.2) implies 
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BOR. = Bn R 1 ?. Since R, C there exists C . e C (R . )n C ( R 1 ? ) . 
1 1 11 i — i—i 


Thus 


R i = {A + BC i | BOR.} = {A + BC i | Bn R^} = R^ 


To verify (3.2) start from 


l^BO £ R^j i d^Bn^ 1 R**) l<j<k-: 


Write P = 


If (3.3) holds with equality for j = l , then 

H M *+1 M 

B n?, R i = Bn 

1 M 

i R ? . Then (3.4) implies 


Bn(P + R i+1 } = BnP + BOR £+l 


so that (Prop. A. 4) 


PD (B + R^ x ) = BnP + POR^ 


Then 


.(3.3) 


(3.4) 


AtPOR^i) c (B + P)H(B + R^ x ) 


c b + Pn R 


M 

l+l 


By Lemma (5.1) of [1] there exists 


M 


C e CfPJnCtR^) 


so that 


R Jl+l = {A + BC I Bn R^JC {A + BC | P} 


C P C N 


l+l 
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in contradiction to (1.5). Therefore (3.3) holds with inequality at 
each j. Since 

k 

d(Bn E RV ) ^ d(B) = k 

i=l 1 

and d(B n ) > 1 (ieJ) the result (3.2) follows. Combining (3.1) 
with Th. 2.1, Corollary, we obtain: 

Theorem 3 . 1 

Let d(B) = k. For the RDP of (1.3) - (1.5) suppose (1.6) is 
.rue, and let = max (I_, (R M )*). There exists a solution 
(E* , , isJ} of EDP if and only if d (E ' ) > n M (l/ M ) , where n M is given 

by (2.5) . 

4. STATE SPACE EXTENSION 
AND POLE ASSIGNMENT 

With the minimal extension of § 2 or § 3 it may happen 
that some poles of the closed loop transfer matrix are necessarily 
fixed at unstable, or otherwise 'bad', locations. It is possible to 
shift the bad poles by additional dynamic compensation. This aim is 
achieved by choosing the extension such that all the fixed eigenvalues 
of A + (B + E ' ) C are 'good'. 

To identify the fixed eigenvalues we need the following Lemmas. 

( emma 4 . 1 

Let 1/e J_, write C = C ( 1/ ) , and let R = max (£, U) . Write 
A c E A + Ce £, and define A c : as follows : if x is the 

coset of x in 1//R , A c x = A c x. Then R and A c are constant with respect 
to Ce C. In particular the characteristic pol y nomial (ch. p.) of 
A c | V has the form ti (X ) tt c (X ) , where tt is the ch. p. of A c and is 
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fixed for all Ce C; tt_ is the ch. p. of A |R; and the roots of tt .can 
— c c c ' c 

be assigned arbitrarily by suitable choice of Ce C. 

Proof 

By [1], Th. 2.3, 

R = {A + BC | BHl/}, Ce C 

and CCC(R). If C,, C, e c and x £ V then A x e V (i = 1,2) and 

“ — it— 

(A -A )x = B(C. - C,)x e BHl/ C R, 

C 1 c 2 1 1 

hence A = A Assignability of the roots of tt follows by [1], 

c i c 2 c 

Th. 2.2. | 

Lemma 4 . 2 

Under the conditions of Lemma 4.1, let a(X) be the minimal 
polynomial of A c# and factor a (A) = oig(A)a b (A), where the polynomials 
a , a b are coprime. Then 

V = R © R g © R fa (4.1) 

where 

R © R g = {x : xel/, a g (A c )x = 0} (4.2) 

and similarly for R © R b * The subspaces R © Rg, R © are fixed with 
respect to Ce C. 

Proof Since a , a. are coprime, 1//R = JT © TT. , where 
— g d g o 


Since R and A c are constant with respect to Ce C, the result follows. 
Lemma 4 . 3 

Let We V and let S = max (C' ,W). Write V £ PW and R £ PS. 

Then 

1. Ve]_, and R = max (C, l/). 

2. V/R = W/S 

3. The fixed eigenvalues of (A + (B + E')C)|W, Ce C' (W) , coincide 
w ith the fixed eigenvalues of (A + BC Q )|l/, C q e £(1/). 

i j roof 

1. If AWCdl + B + E', Al/CPAWCl/ + 8, so l/e U Let R M £ max (C, V ) 
and define 

R° = 0, R u+1 £ 1/ n (AR y + B) (u = 0,1, ... ) (4.3) 

It will be shown that T £ lim R u = R M . Since R y C 1/ and Al/ C 1/ + B, 

AR y C (1/ + 8) n (AR y + B) = R y+1 + 6, 
so that AT c T + 8. Since R y C T C 1/ (y=0,l,...), (4.3) implies 

R W+1 = T H(AR y + B) (U = 0,1,...) 

By [1], Th. 2.1, Tc C_ and TCI/, hence T CR^. On the other hand 

R M = lim R y , where R° = 0 and 

R U+1 = R M n(AR y + B) (y = 0,1,...) 

M A i ta 

Since R Cl/, by induction on y we have R^C R^ , hence R M C T. Thus 

the rule (4.3) computes max (C, l/). 

Applying this result to the pair S,W we have S = lim S y , where 
S =0 and 

S y + 1 = w O (AS y + B + E’) (y = 0,1,...) 
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Thus (Prop. A. 4) PS 


y + ^ = V H (APS y +8) , and comparison with (4.3) yields 
R M = lim R y = lim PS y = PS 
2. In general 

W/S 5 (W + E')/(S + E') © (W n E’ )/(S 0 E') 

But 


W + E ’ (W + E ' ) /E ' 

— . - > -■ ,/c - PW/PS = 1//R 

S + E’ (S + E 1 ) /E ’ 

Also, for Ce C' (W) , 

S= (A+ (B + E')C | (B + E* ) n 0/> 
so that WOE'cS, hence 


(4.4) 


(wnE')/(s n e') * o 

3. Let Ce C' (W) . It will be shown that there exist 

C(S +E')nc(H), C £ £(!/)• (4.5 a,b) 

such that the diagram commutes. By the isomorphisms shown, the result 



W + E '■ 
S + E' 


t 

1//R 


A + (B + E' )C~ 


A + (B + E' )C 


A + BC. 


A + BC, 


► W 

I 


►w/s 


► W + E’ 

s + r 


► U/R 
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will then follow from Lemma 4.1. In the diagram a bar denotes the in- 
duced map in the indicated factor space. Turning to the proof, since 

(A+(B+E')C) SCS, 

• the top square commutes (by definition of bar). Recall that E'niilcS 
and write 

W = ( (S+E')n W) © Z 

= 5©Z (4.6) 

inceAE'= 0 and SeV , there follows S + E'e I/, and there exists 
C^eC (S+E') such that (C^-C)Z = 0. Then A+BC^ is defined, and 

[ A+ ( B+E ' ) C - (A+BC 1 ) ] WCS + E', 

so the middle square commutes. Clearly (A+BC^ (W+E 1 ) C W+E'; since 1/ C 
W+ E\ (A+BC.^1/ = PtA+BC^l/CPdiZ+E') = (/, i.e., C^CdO and (4.5a) is 
true. By (4.4) and (4.6), PZ ■ Z , i.e., 1/ = R©PZ, and C Q exists such 
that 

(C o" C l )R = °' (C o P-C l )Z = 0 (4.7) 

Then 

[ (A+BC q ) P - PfA+BC^JZ = 0 (4.8) 

^ Also, if xeS then PxeR and 

(A+BC q )Px = (A+BC-^Px 

= (A+BC^) (x+e ' ) , for some e'cE' 

= (A+BC^)x+e", for some e"eS+E' 

so that 
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(A+BC Q )Px = P(A+BC q )Px 


= P(A+BC 1 )x + Pe" 


( 4 . 9 ) 


and Pe" eR. Then (4.6), (4.8) and (4.9) imply 


[(A+BC q )P - P (A+BC^) ]1DC ", 


so the bottom square of the diagram commutes. Finally it is clear from 
(4.5a) and (4.7) that (4.5b) is true. § 


We now state a procedure for minimal extension of c.s. R™ to 


achieve both decoupling and an assigned distribution of eigenvalues of 


A+ (B+E ' ) C . We write R^ = R*? and assume the hypotheses of Th. 2.1. 


Extension procedure (EXT) 


Here A will denote the original map in E, not its extension, and 


similarly for B,C. Under the conditions of Th. 2.1, let l^ 1 = max - ( I_, R* ) 
R M 5 max (C, l/ 1 ). For CeC e C(l^), write A = A+BC, and let a(X) 


be the minimal polynomial (mod Rj of A | l/M. Factor a(X) = a (X)a,(X), 

C ^ D 


where the roots of a (a^) are good (bad) . For arbitrary CeC determine 


R M © R = {xrxel^, a (A )xeR M ) 
^ g ' g c 


( 4 . 10 ) 


In (2.4) substitute 1/ = R M @ R^, compute R q = R q (U), and construct a 


minimal solution of EDP as in the proof (direct half) of Th. 2.1. 


With EXT completed, a solution of EDP is now in hand: symbols A 


etc. will again denote the extended maps, defined by (1.8). Write C' = 


0 C' (S. ) , A = A+ (B+E ' ) C . 

1 — l c 


Theorem 4 . 1 


Any solution E', S. (ieJ) of EDP determined by EXT has the followin 


1 . 


(4.11) 


K M C s* C R fl © R g 
2. If CeC ' the ch.p. 7T *^^ of A c l s * can be factored as 


7T* (X) = TT g (X)7T^(X) 


(4.12) 


Here the roots of tt are fixed for CeC* and each root is a root of a ; 

" ^ —————— CJ 

M M 

the roots of tt™ can be assigned as any symmetric set of d(R ) complex 

numbers by suitable choice of C|R M , CeC'. 

3. Write S = S.+...+S,. The ch.p. tt (X) of A I S can be factored as 
1 k — c — c — c 1 


* C (A)= * lc <A) • • • * kC (X) V X) 


(4.13) 


where 


d L - deg TT ic = d ( (R i +R Q ) /R q ) ieJ 


deg tt* = d(R ) 
c o 


(4.14) 


The roots of TT^ c (ieJ) can be assigned as any symmetric set of d^ complex 

M 

numbers by suitable choice of CeC', independent of C|R . 

Proof 

1. By Th. 2.1, EXT determines the S^ such that 

s* = R = n v (R n<R M ©R )) 

1 J 9 

Since N., maximality of the R. implies R . O R^ (jeJ) so that 

! ■ i 

C S* and (4.11) follows. 

2. If CeC', A S*C S*C l^CE, so that A | S* = (A+BC) | S* and C|S* has 

~ c c 

an extension C. : E-+U such that C.eC and A IS* = A |S*. By (4.11) and 
1 l — c^ 1 c 1 

Lemma 4.1 (with R=R M , 1/=^) the ch.p. of A |S* factors as in (4.12), 

C 1 

and the roots of tt™^ are freely assignable by suitable choice of C^| R , 
C l c £' hence by suitable choice of C|S*, CeC’. 
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3. 


The expression (4.13) and assignability of the roots of TT ic follow 

by Th. 1.2 applied to the S^; (4.14) follows by the fact (Th. 2.1) that 

(S.+S*) n E' = 0 and S* = R , hence 
1 o 


(S j .+S*+E')/E' 

(S .+S*)/S* * 

(S*+E ' ) /E ' 


'W/'d' 


I 


Now suppose {E',S^,ieJ} is any solution of EDP, not necessarily 
determined by EXT. Then S*el_ , and PS*eI_. By Lemma 4.3 the fixed eigen- 
values of 


A* = (A+(B+E')C) |S*, CeC ' (S*) , 
coincide with the fixed eigenvalues of 
(A+BC Q )|PS*, C Q eC (PS*) 

As shown in the proof of Th. 1.1, PS^C R^(iR^), hence PS* C R* , and by 
maximality of (/** , PS* c l^. Therefore C Q |PS* has an extension Ac <!/"). 
By Lemma 4.2, 

^ - R"©« g ©R b 

M 

with R 0 Rg given by (4.10). Since the fixed eigenvalues of A* coin- 
cide with those of (A+BC^)|PS*, it follows, if the fixed eigenvalues of 
A* are all good, that 

PS*C V = R M © R (4.16) 

g 

Since the extension constructed by EXT is minimal with respect to the 
properties (2.1) and (4.16), we have proved the following. 

Theorem 4 . 2 

The construction EXT yields a minimal solution of EDP , subject to 
(2.1) and the requirement that the fixed eigenvalues of 
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(A+ (B+E ' ) C ) | S, CeC' 
all be good . 

Remark . Assuming as in [1] that { A | B } = E, we have that {a |8 0E'} = 
E0E'. By the technique used in proving Th. 1.2, it is straightforward 
to show that (E0E')/S can be regarded as a c.s (mod S) for (A,B+E'), 
hence that the only fixed eigenvalues of A are those of A*. 

v C 


c 


5 . EXAMPLE 


. th 


Let n = d(E) =5 and let e^(l<i<5) be the i unit column vector. 


.th 


with 1 in the i row and 0 elsewhere. Let 


A ^ 4 r e^ > e^ r ^ 2 « 6 ^ ] r B 


= row e^, = row e 2 . Writing {•} for the span of the vectors 

bracketed, we have 

^ 2. — ^ ®2 1 ® 3 1 ® 5 ^ * ^2 = ^ ® 2 . r ®3 ' ®4 ' ® 5 ^ 

It is easily checked that 

R 1 = W 2 # R 2 = W l' BnR l = B0R 2 = {e 2 } 


By Th. 5.1 of [1], decoupling by state feedback is not possible. How- 
ever, since ( 1 . 6 ) is satisfied, namely 



R 1 + N 1 - K 2 + W 2 " E ' 


Th. 1.1 asserts that decoupling is possible by use of dynamic compensa 
tion. 

In this example d( 8 ) = 2 = k, and according to §3 any solution 
(i=l,2) of EDP must satisfy 
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PS i = RV i=l , 2 (5.1 

By Th. 3.1 a minimal extension has d(E') = n^l^) given by (2.5), wber< 


l/ M = max (I, (R M )*) 


In this example. 


(R M )* = R^HrM = } 




and one easily computes V - {e^e^. By (2.4), 

R Q (l/ M ) = R^ n R2 n l/ M = l/ M 

Then (2.5) gives n M ( l/ M ) = 1, so that just one integrator is needed to 
achieve decoupling by dynamic compensation. 

To determine the spectrum of A+(B+E')C we follow the procedure EX1 
of §4, and start by finding R M = max (£, l/ M ) . Since B h 1^ = 0, we have 
R M = 0. Since we can take A | l/ M = a| l/*S in our coordinate sys- 

tem 


Al l/ M = 


A 33 A 34 


A 43 A 44 


1 1 
0 0 


so that a(X) = X (X-l) . If unstable eigenvalues are considered 'bad', 

! b 


we have R g = 0 and R K = t/**. Both the bad eigenvalues are fixed in the 


minimal extension determined first. To find the minimal extension 
subject to the constraint that all fixed eigenvalues be good, we set 


(/ = R M 0R g = 0. By (2.4), R q (I/) = 0, and (2.5) gives 


d(E') = n M (l/) = 3 


Exactly three compensating integrators are needed to achieve decoupling 
together with stability of the (extended) closed loop system matrix. 
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The reader may wish to investigate the possibilities with two com- 
pensating integrators. 

6. DECOUPLING AND OPEN LOOP CONTROL 

In previous sections and in [1], the apparently stringent restric- 
tion was imposed that feedback and dynamic compensation be linear. In 
particular the definition of controllability subspace [1] was tied to 
a specific linear feedback structure. We now show that, as regards 
'coupling, nothing is gained by considering more general types of con- 
trol. To this end we show that maximal c.s. can be defined in an open 
loop sense without any assumptions on controller structure. Consider 


x(t) = Ax ( t ) + Bu ( t ) , teT 
x (0) = 0 


( 6 . 1 ) 


on the time interval T = [0,1], and let NCE. Let U denote the class 
of m- vector- valued functions u(.), defined and continuous on T. Denote 
by <*>:TxU->-E the solution of (6.1), i.e., 

<J>(t,u) = f e^ fc S ^ A Bu(s) ds, teT,' ueU 
0 

Theorem 6.1 

Let X be the set of states xeN such that , for some ueU, 

<|>(t,u)eN, teT; <Hl,u) = x 
Then X = R M E max (C,N) . 

U 

Thus R is characterized as the largest set of states in N which 
can be reached from the zero state, by any control whatever, without 
leaving N. 
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Proof Let 


R M = {A+BC | {BK}} 

A A 

and write A = A+BC, B = BK. We claim that 


where 



(R> 



e d-t)A* 


dt 


( 6.2 


(here and below, a prime denotes transpose). In fact zeN(R) implies 
z'e^ 1 fc ^ A B = 0, teT, i.e., z'A^”^B = 0 (j=l,...,n), so ze(R M )^ . Thus 
W(R)C(R M ) 1 , so that R M C {R} , and the reverse inclusion is obvious. 

To show R M CX, let xeR M and note from (6.2) that x = Rw for some 
weE. Set 


v (t) 


B' 


e 


(l-t)A' 

w 


teT 


Then the equation 

A A 

x ( t ) = Ax ( t) + Bv (t) teT 
x (0) = 0 

implies x(T)CR M cW and x(l) = x, where x (T) ={x (t) : teT} . Put 


u ( t) = Kv(t) - Cx (t) , teT 

Then ueU; <t>(t,u)eN, teT; <J>(l,u) = x; and so xeX. 

To show XCR M , let V = max(j[,W) . By [1], Th. 1.1, V « l/ n , where 
V° = W and 

V» +1 = l/ y n a” 1 (l/ y +B) (y=0,l, . . . ) (A^VsfxrAxel/}) 

If xeX then for some ueU (6.1) yields 
x (T) c W, x(l) = x 
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Thus x (T) C V°. If x (T) C V U than x(T) C l/ y , so Ax(T) = (x-Bu) (T)c U U + 6, 
hence 

x(T) C l/ W HA -1 (l/ y +8) = l/ U+1 , 
and by induction x(T)c V. Let CeC(U). Then 
x(t) = (A+BC) x (t ) + Bv(t), teT 
where v(t) = u(t)-Cx(t). Thus 


Bv (T) = (&- (A+BC) x) (T) C V 


so that {Bv (t) } C B fi l/, teT. So, for teT, 


■/: 


x(t) = / exp [ (t-s) (A+BC) ] Bv(s) ds 
0 


e {A+BClBOl/} 


= R 


M 




We now pose an open loop decoupling problem (ODP) as follows. 

Given (6.1), and (1.2) defined for teT, together with arbitrary vectors 
y^eH^(ieJ), find controls u^eU(ieJ) such that 


H ± 4> ( 1 , u ± ) = y ± ieJ 

H^$(T,u^) = 0 i,jeJ; 


(6.3) 


(6.4) 


( ■ 


ier these conditions each u^ affects only the output y^(*)» and y^(l) 


Theorem 6 . 2 

Write ^ = N (H^) (ieJ) . ODP is solvable for arbitrary y^eH^(ieJ) 
if and only if 


R i + W i = E ieJ 


(6.5) 
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where 


( 6 . 6 ) 


R** = max (C, ^ W^) ieJ 

Proof If (6.5) is true than H^R*? = and there is x^eR*? with H^x^ = 

y^. By Th. 6.1 there is u^eU such that 

<^(l,u i ) = x L , <J>(T / u i )c Wj 

i.e. , 

H i 4>(l,u i ) = y L , Hj $ (T ,u^) = 0, j?i 

Conversely if (6.5) fails, then for some ieJ there is yeH ^ such 
M 

that y^H^R?. Therefore any control ueU, such that H^(^(l,u) = y , has 
the property <M1/U)$R*?. By Th. 6.1 

4>(t,u)4 .0. W. 

3 


for some teT; i.e., for this t, Hj<Mt,u) ^ 0 for some jeJ, j^i, and 
(6.4) fails. | 

Comparing Th. 6.2 with Th. 1.1 we have 
Corollary 

ODP is solvable if and only if EDP is solvable, namely if and only 
if (6.5) is true . 

In the definition of ODP the choice U for the class of admissible 
controls, and the choice in (6.2) of common endpoint t=l, are obviously 
not crucial. In fact we have shown implicitly that a wide class of 
dynamic decoupling problems is equivalent to the EDP of §1. 
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CONCLUDING REMARK 


Taken with its predecessor [1], the present article provides 
effective machinery for the formulation and solution of the de- 
coupling problem. The results prescribe the synthesis of dynamic 
compensation by which decoupling can be realized, and clarify the 
conditions under which such compensation exists. Nevertheless, 
further aspects of the problem remain for investigation. These 
include computer implementation, sensitivity analysis, and perhaps 
most important, a deeper account of algebraic structure. 


APPENDIX 


We collect here some auxiliary results; verifications, when 
straightforward, are omitted. 

1. Let l/^(ieJ) be arbitrary subspaces. Let 

v! = £ v., v' I n V* 

1 if>i 5 i 1 

Then 

v* =£ i'. n v* =2 v. n i/* 
i 1 i 1 1 

- E i'!/ jcJ 

1 1 

2. If u i - + v * (ieJ) then u * = l/ * 

3. if X = n 2 l/.n / for some V , then 

i j^i 3 

x=nE u.nx (a. 3) 

i j*i 3 

4. By definition the l/^(ieJ) are mutually independent if and only if 

* * 

V = 0, i. e. , O (/^ » 0 , ieJ. More generally: 

Prop, A. 1 

V* is the smallest subspace 1/ such that the factor spaces (l/^+ 
l/ Q )/l/ o are independent . 

Proof Independence of the factor spaces is equivalent to 

V Q =E(V t 'o ) n(t 'i +t ' 0 ) <A.4) 


(A. 1} 
(A. 2) 
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* u 

From (A. 4), V = lim 1/ (y=0 , 1, . . . ) , where 
l/° = 0, l/ w+1 (l/. + U M ) n (l/*+l/ u ) 


(A. 5) 


By (A. 2), V satisfies (A. 4), and (A. 5) implies that any solution \J q of 
(A. 4) contains 1/ . | 

5. By Prop. A.l, 


d ( l/ j ,/ ( l/ i 0 l/*) ) =53d((l/. + l/*)/l/*) 
i i 

= d(53 (l^ + l /*)//) 

= d(J3v.)/V*) 


so that 


A [U i ,cr ] =53d(i/ i ) - d(Ei/ i ) = a < t/ ± n v*) - d(i/*) (a. 6) 


6. If U , V , W are arbitrary subspaces, 

( u+w ) n (t/+g ;) ^ (ti+t/) n u 

• uni/ + w unw + i/nw 


( 


(A.7) 


Prop. A. 2 

Let S . ( i e J ) , 1/ , E ' be such that l/H E' = 0, S C 1/ © E ' . Define 


S e£s 
i 


l 

and 


s e n 2-f (S .+e 1 ) n (1©E') 

° i f 7 i 3 


Then 
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d(E') - 6 ^ + & 2 + P 

where 

6 1 = M(S..+S 0 +E')/(S 0 +E'), J] 

6 2 = A[ (S.+E') n (S o +E')+S*)/(S*+E') , J] 

P =Ed[((S.+s*) n E')/(S. n E'+s* n e') ] 
i A 

+ 2d[ (s. n E')/(s . n s* n e' ) j 

i x 

+ d (S* fl E') + d (E '/ (S HE' ) ) 


(A. 8) 

(A. 9) 
(A. 10) 

(A. 11) 


Proof . The proof is a direct computation, starting from the easy 
identity 

d(E') = d (S) -d ((S+E'J/E*) + d(E'/(S n E') ) 
and using (A.l) - (A. 7); from (A. 3) note especially 

s n = n 2 (S.+E*) n (S +E') ■ (A. 12) 

° i j^i D ° ■ 

8. Prop. A. 3 

If in (A. 8) , d(E') = 6^ then S* + E' = S Q , (Sh-S*) O E' = 0 (ieJ) , 
and E'c S. 

* 

Proof , p = 0 implies SHE' = E ' , i . e . , E'c S; also S HE' = 0 , hence 

* 

S^PE' = 0 (ieJ) ; so that, from the first summation in (A. 11), (S^+S )OE 
= 0 (ieJ) . Also, 6 2 =_ 0 implies that the bracketed factor spaces 
in (A. 10) are independent; by Prop. A.l and (A. 12), 

s* + E'd n 2 ( (S.+E ' ) n (s +E')) = s n (a. 13 ) 

i j^i 3 ° ° 
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By (A. 3) and the definitions of S , S q , 

s d n E s . n s* = s* , 
i j^i 3 

* 

hence the reverse inclusion holds in (A. 13), so S + E' = S Q . 
9. Prop. A. 4 

For arbitrary U, \J , W, if 1 

u n (i/+w) = ar. u + unw 

then 

\j n (u+w) = uni/ + i/nw 

10 . Prop. A. 5 

For arbitrary U, V and a map T, 

T (Uni/) = (TU) O (Tl/) 

if and only if 

(u+i/)nw(T) = unw(T) + i/nwcn 
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